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Abstract 

Kamenev and Mezard, and Yurkevich and Lerner, have recently shown 
how to reproduce the large-frequency asymptotics of the energy level 
correlations for disordered electron systems, by doing perturbation the- 
ory around the saddles of the compact nonlinear a model derived from 
fermionic replicas. We present a critical review of their procedure and 
argue that its validity is limited to the perturbative regime of large fre- 
quency. The miraculous exactness of the saddle-point answer for /3 = 2 
(unitary symmetry) in the universal limit, is shown to be a special feature 
due to the Duistermaat-Heckman theorem. 

1 Introduction 

Suppose we are given a disordered statistical mechanical system with free energy 
/ = — InZ, and we are to compute the disorder average (/). To tackle that 
task, it is popular among theoretical physicists to follow a recipe called the 
replica trick. Instead of attempting to calculate (InZ) directly, one computes 
the disorder average /„ = (Z n ) for all positive integers n = 1,2,.. .00, which 
is be done by introducing n copies, or replicas, of the system. Then, recalling 
the identity \nZ — ]im n — > o(Z n — l)A l i an d assuming some extension of the 
discrete set (/ n ) n eN to an analytic function f(u) on C, one hopes that /'(0) 
equals — (In Z) . An obvious problem with this trick is the non- uniqueness of 
the analytic continuation. Indeed, if f(u) satisfies f(u)\ — fn> then so 
does f(u) + c/r(— u), where T(u) is the Euler gamma function and c is an 
arbitrary constant. Notice that the increment c/T(— u) vanishes for all positive 
integers u — n but has finite derivative at w = 0. Thus we do not know 
whether the answer for — (lnZ) is /'(0) or /'(0) — c. To determine the unknown 
constant, which may be a function of temperature, disorder strength etc., one 
must have some control on (Z u ) , for example a bound on its asymptotic behavior 
as u — > ioo. Sadly, such control is not always available. 
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Some time ago, Verbaarschot and the author (VZ) addressed JO the above 
concern in the context of disordered electron systems, i.e., the problem of com- 
puting averages of retarded and advanced Green's functions G ± (E) = (E ± 
ie — H)^ 1 of a disordered single-particle Hamiltonian H. For such systems, one 
implements the replica trick by starting from a generating function 

Z m , n (oj, Ji, J a ) = (Bct m (E + uj + ie-H + Jx) Det n (E - ie - H + J 2 )) 

depending on matrix sources J\ , J 2 . Angular brackets denote averaging over 
disorder. For positive integers m and n, this generating function can be ex- 
pressed as a Gaussian integral over anticommuting complex fields ("fcrmionic 
replicas"), and for negative integers over commuting fields ("bosonic replicas"). 
VZ focused on the spectral two-point function, 

S 2 {uj) = (tt{E + uj + ie - H)- 1 Ti{E - ie - H)- 1 ^ (e > 0), 

for the simple case of Wigner-Dyson statistics for j3 — 2, i.e., the Gaussian 
ensemble of complex Hermitian matrices with unitary symmetry (GUE). By 
following the standard procedure @, ||, of mapping the problem on a nonlinear 
a model, and then doing the natural analytic continuation m = n — ► 0, VZ found 
that fermionic and bosonic replicas give different answers, and both answers 
differ from the exact result. What is reproduced correctly by the replica trick is 
the asymptotics of ^(w) as lo — > 00. This finding led to a general consensus that 
the replica trick for disordered electron systems is limited to those regions of 
parameter space where the nonlinear a model can be evaluated perturbatively. 

This consensus is now being challenged by recent proposals of Kamenev and 
Mezard d, ^| , and of Yurkevich and Lerner . A partial list of their claims 
is as follows: i) According to Kamenev and Mezard ||, past work using the 
replica trick for disordered electron systems suffers from the tacit assumption of 
the absence of "spontaneous breaking of replica symmetry" . ii) When replica- 
symmetry broken saddle points are included, the "standard" nonlinear a model 
derived from fermionic replicas gives the correct oscillatory asymptotics of the 
level correlation functions, iii) Perturbation theory in uj~ l yields a "systematic" 
expansion for the amplitudes of the oscillatory terms, iv) The fermionic replica 
trick with replica-symmetry breaking is "nonperturbative" and reproduces the 
results of the supersymmetric method, v) Yurkevich and Lerner speculate that 
if limits are taken in the proper order (first n — > 0, then ui — > 0), the fermionic 
replica trick may reproduce also the singular part of the two-point function. 

These claims and speculations add up to a substantial revision of the ac- 
cepted picture. They provoke some sort of response, and it is the purpose of 
the present paper to offer clarification. Section ^ calculates in some detail the 
one-point function, Section ^ reviews the two-point function, Section |4| explains 
why replicas can be manipulated to yield the exact answer for the spectral 
correlations of the GUE, and Section summarizes the conclusions. 
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2 Fermionic replicas: one-point function 



Let us start from the beginning and consider first the spectral one-point function 
Sx(z) — (Tr {z — H)^ 1 ) . The fundamental difficulty plaguing fermionic replicas 
for disordered electron systems is already visible in this case. 

On elementary grounds, the function S\(z) has a branch cut along the real 
axis (or the line of spectral support of the Hamiltonian H ) , and its discontinuity 
across the cut yields the average density of states, p(E): 

p(E) = — lim (SAE- ie) - SUE + ie)) = - lim Im SUE - ie) . 

The replica trick, instead of attempting to calculate Si directly, looks at the 
discrete family of functions 

f n (z) = (Tr (z - H^Be^iz - H)) (z £ C) 

for n £ Z. The member fo(z) of this family is the desired quantity S\(z). In 
the fermionic version of the replica trick, one computes f n {z) for all positive 
integers n = 1, 2, . . . oo, and hopes to infer S±(z) by setting n = at the end of 
the calculation. On brief reflection, however, this hope must appear ill-founded. 
Indeed, it is easy to see that, if all moments {(Tr i? p )} p= i ) . --in _i exist, the 
function f n {z) is analytic in z for all positive integers n £ N, which means in 
particular that f n (z) is continuous across the real axis, and 

Pn {E) = lim (f n (E - ie) - f n {E + ie)) = - lim lmf n (E - ie) = 

is zero, for all E £ K and n £ N. This is beginning to look bizarre. Given 
that all imaginary parts p n {E) vanish identically for n = 1,2, ...oo, how can 
we hope to predict, based on these data only, the nonzero value of the average 
density of states p(E) = p n (E)\ n _ ?! 

Let us postpone the answer for a moment and jump to the computation of 
fn{z) (n > 1) using fermionic replicas. For simplicity we focus on the case of 
N x N random matrices with (3 — 2 (unitary symmetry) and take the probabil- 
ity measure for the complex Hermitian matrix H to be Gaussian with weight 
function P(H) = const x exp(— NTrH 2 /2w 2 ). By a standard sequence of trans- 
formations H [jj, the function f n (z) can then be cast in the form of an integral 
over u(n), the Lie algebra of the unitary group U(n): 

/„(*) = - / Tr(z-g)- 1 Dct Ar (z-Q)e JVTr( 5 2 / 2 - 2 dg, 
n Ju(n) 

where dQ is the flat measure on u(n) normalized by f u ( n } eNTl ® ^ 2w dQ = 1- 
The elements of u(n) are skew-Hermitian matrices, whence Tr Q 2 = — Tr Q^Q < 
and the integral exists. Note that the integrand is invariant under the adjoint 
action Q i— > UQU^ 1 of U(n) on its Lie algebra. 
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For large N, the integral over Q is saturated by the solutions of the saddle 
point equation 

— (TrQ 2 /2«; 2 + Trln(z - Q)) = = Q/w 2 + (Q - z)" 1 . 

This equation is of real type. Its solutions therefore come in Hermitian conjugate 
pairs: if Q is a solution for z — E £ M, then so is Q\ In particular, for \E\ < 2w 
there exists a solution 

Q = we ie xl ni e lS = E/2w + iy/l - (E/2w) 2 , 

which is proportional to the unit matrix, and is paired with = wc _lfl x 1„. 
These two solutions are unique in that they are stable under the adjoint action 
of U(n). For n > 2, there exist also solutions with less symmetry, which organize 
into smooth saddle-point manifolds. They can be constructed by starting from 
a diagonal solution with n — p (resp. p) matrix elements equal to we l6 (resp. 
we~ l6 ), and then applying the adjoint action of the symmetry group U(n). 

The point to be made here is this: proper application of the saddle-point 
method to the integral representation for f n (z) requires us to sum over pairs 
of conjugate saddles. It is clear that the contribution from a saddle point (or 
saddle-point manifold) Q is precisely the complex conjugate of the contribution 
from its Hermitian conjugate Q\ if z = E± iO. Therefore, we can be sure that 
the result for f n (z) on the real energy axis is real: 

/„(E) = ^(E) + ^(£)el, 

and p n {E) = ir^ 1 lim £ ^o f n (E — ie) is zero for all n = 1, 2, ...oo. This, of course, 
was to be expected in view of what we saw at the beginning of the section. 

We now return to the question posed before. Since p n {E) vanishes iden- 
tically for all n £ N, how can we manage to produce a reasonable result for 
p(E) — p n (E)\ _ Q ? Clearly, the answer is that some additional input must 
be injected. What is missing from the formalism with fermionic replicas is the 
entire information about causality, i.e., the distinction between retarded and 
advanced Green's functions (z = E + ie, versus z = E — ie). This information is 
present in f n (z) only for n < 0. In the bosonic replica trick, where one evaluates 
f n (z) for n = —1, —2, . . . — oo, causality is reflected in the Q-integral representa- 
tion as a singularity of the integrand, which restricts the number of saddles that 
are accessible by continuous deformation of the integration contour Js) . It turns 
out that for all n < 0, there exists just one saddle point which is accessible, and 
all others are inaccessible. On the other hand, in the rigorous formulation of 
fo{z) by the supersymmetric method Q, there are four saddle points, two acces- 
sible and two inaccessible ones. Both the supersymmetric formula and bosonic 
replicas will be discussed in more detail below. 

2.1 The recipe 

As we have seen, the reason for the vanishing of p n {E) [n £ N) is the pairwise 
appearance of saddles as conjugate pairs It is therefore clear that, 
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in order to assist the fermionic replica trick and extrapolate p n (E) to a non- 
vanishing answer for n — 0, we must come up with a recipe for breaking the 
symmetry between the members of conjugate pairs. Let us denote by Q p the 
diagonal matrix with the first n — p entries equal to we and the last p entries 
equal to roe -19 . If z = E G R, the matrix Q P lies on the saddle-point manifold 
UQn-pU^ 1 (U G U(n)), so the pairing of saddles is p <-> n — p. 

For definiteness, let us fix z = E — iO, in which case ImSi(z) > 0. According 
to the saddle-point equation (z — Q)^ 1 = Q/w 2 , the expression 

lm(E -we 10 )- 1 = Imc t8 /w = V 1 = (E/2w) 2 

is positive for \E\ < 2w, which is known [1(J to be the range of the energy 
spectrum in the limit N — > oo. Since the integral representation for f n (z) 
contains the factor Tr (z — Q)" 1 , the "good" saddle point, contributing to f n {z) 
with an imaginary part of the desired positive sign, is Qo = we' 8 x l n . By the 
same token, Q n = we~ l9 x 1„ is a "bad" saddle point. Alternatively, we could 
say that Qo is "causal", whereas Q n is "acausal". By heritage, or continuity, 
the natural recipe now is to declare the saddle-point manifolds UQpU^ 1 with 
p = 0, 1, 2, ... to be "causal", and their conjugates p = n,n — l,n — 2, . . . to 
be "acausal" . (Clearly, this only makes sense if there is some mechanism that 
terminates the sequences before they overlap.) Anticipating the emergence of 
causality for n — 0, we then retain only the contributions from the "causal" 
saddles, while throwing away the others. The resulting fake expressions for 
f n (E — it) have nonzero imaginary parts, and extrapolation to n = can now 
give a finite density of states. The result so obtained for \E\ < 2w reads || 

p = {N/irw) sin6 + O(N ) + ... 

-(Attw sin 2 9)' 1 cos (n(29 - sin 20)) +C(A" 1 ) + ... . (1) 

The leading term is Wigner's famous semicircle law, sin# = y/l — (E/2w) 2 , and 
stems from the invariant saddle point p = 0. The subleading (oscillatory) term 
originates from the saddle-point manifold p = 1, which is generated by applying 
the adjoint action of U(n) to Qi, the diagonal matrix with n — 1 entries we 16 
and one entry we~ l9 . The stability group of Qi in U(n) is U(n — 1) x U(l). 
Hence the saddle-point manifold p = 1 is isomorphic to the quotient U(n)/U(n— 
l)xU(l). More generally, the saddle-point manifold with index p is isomorphic 
to the Grassmannian G UtP = U(n)/U(n — p) x U(p), which is a symmetric 
space of type Alll in Cartan's classification . Carrying out the saddle-point 
approximation requires integrating over the directions transverse to the saddle- 
point manifold approximately, and over the saddle-point manifold itself exactly. 
The latter integral produces a volume factor vol(G„. p ). For p = 1, this volume 
goes essentially as r(n)^ 1 (see for example the appendix of fj|), which combines 
with the prefactor n" 1 of the integral representation for f n (z) to produce a finite 
limit n _1 r(n) _1 — > 1 for n = 0. More generally, "causal" analytic continuation 
[^| of the volume of G„ jP yields a power n p near n = 0. The contributions from 
the saddle-point manifolds p > 2 therefore extrapolate to zero at n = 0. 
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Remarkably, the above TV -1 expansion of the average density of states agrees 
[^| with the asymptotic behavior that can be inferred from an exact representa- 
tion in terms of orthogonal polynomials. In view of the rather ad hoc nature of 
the derivation, mathematicians will throw up their hands in horror or despair, 
while physicists are much intrigued. In any case, putting aside questions of 
mathematical rigor, it is definitely desirable to gain a better understanding of 
why the procedure works. 

2.2 Supersymmetry 

One avenue towards some understanding is to make a comparison with the 
rigorous representation |], ^| of Si(z) as an integral over a 2 x 2 supermatrix 



Here DQ stands for the flat Berezin form DQ = (2ir)~ 1 dxdy d 2 /d£dr], and the 
integral is over x S K and iy £ iR. The symbols STr and SDet denote the 
supertrace and the superdeterminant. 

The above Berezin integral for Si(z) has a total of four saddle points, ob- 
tained by going through the various sign choices in Q — w diag (e ±l9 ,e ±ie ). Two 
of these are "bad" , as they cannot be reached by continuous deformation of the 
integration contour for x without crossing the iV-th order pole of $DeXT N {z-Q) 
at z = x. As before, let us fix z = E — iO. Then the accessible saddle points 
are Qo = wdiag(e 10 ,e i9 ) and Q\ = w diag(e* e , e~ l9 ). The first one gives rise Q 
to Wigner's semicircle, and on comparing with the fermionic replica trick, one 
begins to suspect that perturbation theory around the second one might yield 
the oscillatory contribution to the asymptotic result (^). I have checked that 
this is precisely the case. 

In hindsight, the coincidence is not all that surprising. At the non-invariant 
saddle point Qi = w diag(e l6) , e~ 10 ), both £ and 77 are Goldstone fermions of the 
supersymmetric integrand SDet _Ar (z — Q) exp — A^STr Q 2 /2w 2 . On the other 
hand, the p = 1 saddle-point manifold U(n)/U(n — 1) X U(l) of the fermionic 
replica trick is isomorphic to the complex projective space CP n _i, which is a 
manifold of real dimension 2{n — 1). At n — 0, this has real dimension —2, and 
has to be understood as a space generated by two anticommuting degrees of 
freedom, say £ and 77. Hence, integrating over CP n _i ~ U(n)/U(n — 1) x U(l) 
and then setting n — is just a complicated way of taking two derivatives 
d 2 /d£,dr). The latter is what the Berezin integral formula does. The agreement 
between the asymptotic expansions to leading order motivates us to propose 
the following conjecture: the fermionic replica trick, augmented by the plausible 
recipe of retaining only the saddles p — and p — 1, is perturbatively equivalent 
to the supersymmetric integral representation for S% (z). In other words, pushing 
the asymptotic expansion around the two saddles Qo and Qi to higher order, 
we expect to get agreement in every order of iV _1 . 




Si(z) = — DQTr(z-Q)- 1 ST)ct- N {z-Q)e 




-NSTiQ 2 /2w 2 
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Can this conjecture be put on a more solid footing? With the goal of justi- 
fying better the asymptotic expansion obtained from the fermionic replica trick, 
we are now going to take a more close-up look at the procedure of analytic 
continuation in n. 

2.3 Analytic continuation? 

To keep the calculations as simple as possible, we start from the generating 
function 

F n (z) = (pet n {z-H)) , 

aiming to eventually apply n~ 1 d/dz and analytically continue to n = 0, so as to 
extract Si(z). As we have seen, the values of F n (z) at the positive integers n do 
not contain sufficient information to permit analytic continuation. Moreover, it 
is not known how to compute F n (z) for n Z, at not least not directly. We 
do have access, however, to F n (z) at the negative integers n, by means of the 
bosonic replica trick. Let us see what additional information we can collect from 
there. 

By some elementary transformations, we arrive at the identity 

/ Dct- N (z-Q)e- NTl Q2 / 2w2 dQ (n e N) . 

J iu(n) 

We are now integrating over the space of Hermitian matrices, viewed as iu(n). 
The saddle-point equation governing the large- N limit still reads [z — Q)^ 1 — 
Q/w 2 and, as before, the invariant p = saddle point yields [l4| Wigncr's 
semicircle. Other saddle-point manifolds with index p = 1,2, ...,n exist, too. 
The difference from before is that these are now disconnected from the integra- 
tion domain iu(n) by the singularities of the function Det (z — Q). In view 
of what happend for fermionic replicas, one might think that one should ignore 
disconnectedness and incorporate the other saddles. Surely, moving across the 
iV-th order pole of T)et~ N (z — Q) just produces some residue, and one could 
argue that this is schematically given by S N_1 /SQ 1 ^ -1 exp(— iVTr Q 2 /2w 2 ) \q_ z , 
which is exponentially small and hence negligible, for large N. However, we are 
on the wrong track. A little thought shows that the saddles p = 0, 1, ...,n are 
aligned transversely to the contour of integration for Q. This makes it impossi- 
ble to deform the contour so as to pass through more than one member of the 
sequence of saddles without backtracking. (In contrast, for fermionic replicas 
the alignment is longitudinal, and a constant shift Q — > Q + E/2 puts all the sad- 
dles right on the contour.) Moreover, if one proceeds anyway and computes the 
contribution from p — 1 , one obtains an expression containing the exponential 
factor 

exp (nN(i6 - e lt> /2) - iN{29 - sin 29)) . 

(Again, this is for z = E — iO.) Near n = this can be seen to increase without 
bound as z moves away from the real axis toward z — —ioo. Such behavior is 
incompatible with analyticity of the Green's function at infinity. For all these 
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reasons, it is safe to say that F n (z) at the negative integers n is exhausted by a 
single saddle point, p = 0. Evaluation of the corresponding integral in Gaussian 
approximation gives 

F n (2wco S e - iQ) = w nN ggir gjg ~ l 3 + 0(N -i) (_ n€N ). ( 2 ) 

Vl-e 2ie " 

This result will place a strong constraint on the analytic continuation of F n . 

Next we return to fermionic replicas, and carry out a more complete saddle- 
point evaluation of 

F n {z)=l Vet N (z~Q)e N ^Q 2 / 2w2 dQ 
Ju(n) 

for n G N. As before, let Q p denote a diagonal matrix with n — p entries 
we 1 ® and p entries we -1 *. By acting with the symmetry group U(n) on Q p 
we get the saddle-point manifold G n . p = U(n)/U(n — p) x U(p). To perform 
the saddle-point approximation, we need to factor the integration domain near 
Gn tP in a suitable manner. For this purpose, let u(n) be decomposed into two 
linear subspaces: u(n) = V + W, where the elements of V = u(n — p) + u(p) 
commute, while those of W anticommute, with Q p . We then introduce adapted 
coordinates by 

Q = U(Q P + X)U- 1 , 

where UQpU^ 1 parametrizes G ntP , and X € V spans the directions transverse 
to G„ iP . By an elementary calculation, the volume element dQ under this 
factorization transforms into 

Det &d{Q p + X)\ w dX dU . 

Here, dU is a U(n)-invariant measure for G„ jP , dX is a flat measure on V, 
and ad(y)| w denotes the commutator (or adjoint) action of Y € V + iV on the 
elements of W + iW. By expanding the integrand T>et N (z — Q) exp TVTr Q 2 /2w 2 
w.r.t. X and doing the integral over X in Gaussian approximation, we obtain 



F n ~ ^Det-^Qp/w 2 ) exp(ATrQj;/2«; 2 ) Detad(Q p )| 

p=0 



if 



X 



u;" Jv ^vol(G„, p )(7V/2^)^-rt 1 ([1 p) / 



P =o VI - e 2 ^ v w Vl-e- 2 '^ 

x exp (nN{\e 2ie - iff) + ipN(29 - sin 26»)) (n € N) . (3) 

An analogous expression for the derivative n~ 1 dF n /dz was obtained in [EJ. We 
now face the task of finding an analytic continuation of F n to n € C which 
respects the results (||) and (^) for n G Z. Kamenev and Mezard j5j have 
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suggested that the solution to this problem is to analytically continue vol(G„ iP ) 
as 

vol (g„, p ) = e^-*) n m , 

and extend the upper limit on the sum over p in (^) to infinity. The latter does 
not change the value of F n for n € N, as T(n — p + vanishes for p > n + 1. 
For the same reason, the terms p > 1 vanish at the negative integers n, in 
agreement with the constraint posed by (|2|). So far so good. Unfortunately, 
when n is not an integer, the infinite sum over p is ill-behaved. The culprit is 
the rapid growth with p of the (analytically continued) volume factors: 

{2w) 2 P -n-i = T ® = T(p)T(p n) ^ S mn(p - n) , 

vol(G n , p _i) T{n-p + l) 

which makes the sum over p diverge. There appears to exist no obvious remedy 
for this difficulty, and we have to concede that extension of the sum over p to 
infinity does not define F n {z) away from the integers. Thus we have so far failed 
in our attempt to perform a bona fide analytic continuation, and the recipe of 



Section 2.1 still awaits justification. 

Nevertheless, there seems to be some truth to the proposal by Kamenev and 
Mezard. If we boldly ignore the divergence issue, and formally apply {^Kn)~ x djdz 
at n = 0, then all the terms for p > 2 disappear, and the answer from p = and 
p = 1 can be shown to agree with the large- N asymptotics extracted from the 
exact formula for S±(z) as a Bcrczin integral. Upon further thought, the failure 
of the proposed analytic continuation can be interpreted as follows. Recall that 
the failure is caused by the late terms in the sum over p, which in turn stem 
from our insisting that the analytic continuation reproduce the result ([|) for 
all n £ N. But at fixed N < oo, the saddle-point approximation leading to 
that expression breaks down when n becomes large. Indeed, for n 3> N the 
Q-integral formula for F n (z) is inappropriate, and it is better to return to the 
original representation 

F n {z) = J Det"(z - H) e~ NTl - H2 / 2w2 dH . 

By evaluating this integral in saddle-point approximation, one gets a an answer 
formally similar to (j^) . A major difference is that the sum over saddle-point 
manifolds p now terminates at N, not n. Thus, making n bigger no longer 
extends the range of p. The lesson to be learnt from this is that the approxi- 
mate result (|3|) must not be trusted for large values of n. It is likely that the 
divergence of the sum over p is fake and is somehow cut off when a more correct 
approximation for large n is used. When N is big enough, this cutoff will have 
a negligible effect on the extrapolation to n — 0. On the other hand, for small 
values of N the cutoff will probably lead to some unknown correction. 

We could have anticipated the existence of such corrections already from 



our simplified presentation in Section 2.1. Surely, the recipe of throwing away 
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half of the saddles makes sense only in the perturbative domain of isolated 
saddles. When the saddles are not well separated (as is the case for small AT), 
their asymptotic series in some sense "interfere" , which is likely to give rise to 
corrections of exponentially small type ~ e -c°nstxiV_ p or ^he one-point function, 
the existence of such corrections is of minor concern, as we are usually satisfied 
with knowing the large- TV asymptotics. As we shall see, however, the situation 
is less favorable for the two-point function. Note also that the question of what 
happens beyond the saddle-point approximation or, more precisely, whether a 
non-analytic function can be reconstructed from its asymptotic series (by Borel 
resummation, for example), is a very difficult issue. 

Let me finish the section with a word on vocabulary. In my opinion it is inap- 
propriate to call the appearance of the p — 1 saddle-point manifold an instance 
of "spontaneous breaking of replica symmetry" . Spontaneous symmetry break- 
ing, a phenomenon which may occur in physical systems with infinitely many 
degrees of freedom, means that a symmetry of the Hamiltonian is not manifest 
in the (ground) state of the system. The symmetry group of the present problem 
is U(n). The saddle-point manifold CP„_i ~ U(n)/U(n— 1) x U(l) is a beauti- 
ful symmetric (and replica-symmetric) space, whose geometry is invariant under 
the action of U(n). It is true, of course, that an individual point on the manifold 
is not invariant under U(n). However, our job is to integrate, and therefore the 
issue is not invariance of individual points but invariance of the domain of inte- 
gration. The integration manifold CP„_i ~ U(n) /U(n— 1) xU(l) definitely does 
not break U(n) replica symmetry, just as the familiar two-sphere in Euclidean 
3-space, S 2 = SO(3)/SO(2), does not break SO(3) rotational symmetry. 



3 Fermionic replicas: two-point function 

Let G be one of the compact matrix groups Sp(2m + 2n), U(m + n), or 0(2m + 
2n). Acting with G on a fixed matrix £3 by conjugation, we get an orbit 
of elements we denote by Q = U^U^ 1 (U e G). For G — Sp, O we take 
£3 = diag(l 2m , -l2n); and for G = U, £3 = diag(l m , -1„). To keep the 
dimensions explicit, we write X m n for the orbit of G on £3. The stability group 
of £3 is a subgroup, H , of G. Letting G act on £3 by conjugation amounts to 
taking the quotient of G by H, so AT TOj „ is isomorphic to the coset space G/H 
(which is moreover a symmetric space). 

Now let dQ be a G-invariant measure on X m n , and consider the generating 
function 

Z m ,n(") = J exp^Tr £ 3 <2) cZQ. 

This is the zero-dimensional limit of the nonlinear a model of Wegner and 
Efetov Q, truncated to the compact sector arising from replicated fermionic 
fields ||. The generating function for G = Sp, U, O is intended to describe the 
universal low-frequency limit of a disordered electron system with symmetry 
index f3 = 1,2,4 (i.e., orthogonal, unitary, or symplectic symmetry), in this 
order. The dimensionless parameter u is a frequency measured in units of the 
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mean level spacing divided by n, and a = 2, 1, 1 for /3 = 1, 2, 4. For future use, 
note that the integrand is invariant under the action of the little group H. 

The transformation Q — Q is a natural isomorphism from X m ^ n to X n ^ m . 
Using it in the expression for Z m ,„(w) we obtain the symmetry relation 

Z m ,n(^) = Z n ^ n ( — L0). 

This symmetry does not come as a surprise but reflects the fact that Z m , n (uj) 
is a low-frequency approximation to 

(Dct m (E - ie + u/2 - H) Det"(£ + ie - uj/2 - H)} . 

Because this is an analytic function of uj/2 — ie for m, n e N, nothing is lost by 
setting e = 0. We then see at a glance that the generating function is invariant 
under the combined operation of exchanging m <-» n and sending w — > —ui. 
Recall the definition of the two-point function 

S 2 (u,v) = (tt (u - H)- 1 Tr (v - H)- lX J . 

We are interested in the case where the arguments of S2 he on opposite sides of 
the real energy axis: u = E\ — iO, v = E 2 + iO, and their difference is a fixed 
multiple oj/tt of the mean level spacing. The resulting function is still denoted 
by S 2 {u>). 

In the fermionic replica trick for the two-point function, one hopes to extract 
5*2 (ui) from the knowledge of Z m . n (uj) for all positive integers m and n, by 
applying two derivatives — (mn)~ 1 d 2 /dco 2 and setting m = n = 0. The bad news 
is that this hope is frustrated by the symmetry relation Z min (u>) = Z ntm (—u>). 
If the replica limit exists in a naive sense, then it cannot matter whether we first 
set m = and then n — 0, or the other way around. Another option is to set 
m — n = first and m + n = afterwards. If we proceed in the latter fashion, 
the symmetry of the generating function reduces to Z n>n {u) = Z n _ n (—uj). Thus, 
Z n ^ n {ui) is an even and hence real function of u>. By naive extrapolation to n = 0, 
we would conclude that S^w) is an even and real function of its argument. But 
this clearly is nonsense, for S2 is neither even nor real. For example, if Imu) < 
the universal answer for the case f3 = 2 is S^w) = 1 — 2ius~ 2 c~ tu sinaj. 

The reason for the nonsensical answer produced by our use of fermionic repli- 
cas is the same as for the one-point function: the analyticity in ui of the func- 
tion Z m ^ n (iS) means that the fermionic replica trick suffers from the deficiency 
of being entirely ignorant of causality, or the distinction between retarded and 
advanced Green's functions. To get a reasonable answer for £2(0;) and other 
quantities, we need some recipe for adding the missing information. Such a 
recipe is available in the pcrturbative domain of large u), where the integral for 
Zm,nW) can be evaluated by stationary-phase approximation around the saddle 
points of the integrand. As in the case of the one-point function, the recipe is to 
select those saddle points that contribute with the behavior dictated by causal- 
ity. These are the "good" or "causal" saddle points, and they are retained. The 
remaining ones are "bad" or "acausal" and are thrown away. 
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Let us review in somewhat more detail how this works. For concreteness, 
we fix the sign of the imaginary part of w to be negative. The saddle-point 
equation for the integrand exp(iwTr T, 3 Q/2a) is easily seen to be [Q,S 3 ] = 0. 
This equation is solved by the matrix Q = S3 in all cases (3 = 1, 2,4. For the 
sign choice made, this is a good saddle point, as is suggested by the observation 
that the saddle-point value e iajTrS 3 is a maximum of the integrand for Imw < 0. 
(Actually, the true reason is that the saddle-point value e ll ^( m + n ) r e 2l ^( m + n ) 
decreases in magnitude when lu is moved into the lower half of the complex plane, 
for negative m + n. Growth would be incompatible with the analytic properties 
of the Green's function.) By the same token, the saddle point Q — —S3 is bad 
and is thrown away. When the sign of Imw is changed, the roles of Q = £3 
and Q = —S3 are reversed. Note that the two matrices Q — ±£3 are special 
in that they are invariant under the action of the symmetry group H , which 
means they really are isolated saddle points. 

Other solutions Qi, Q2, ■■■ of the saddle-point equation are obtained by 
starting from the good Qq = S3 and exchanging the positions of a small number 
p (or 2p) of the m (or 2m) entries +1 with the same portion of the n (or 2n) 
entries —1. Since they descend from a causal parent, the saddles so obtained are 
expected to be still "causal" and are retained. Their negatives —Qi, —Q2, ■■■ are 
thrown away. Unlike Q a , the diagonal matrices Qi,Q2, ... are not stable under 
the action of the symmetry group H but belong to saddle-point manifolds or 
orbits generated by H . Once again, these manifolds belong to the category of 
symmetric spaces. For example, the orbit of H = U(m) x U(n) on Q\ for (3 — 2 
is 

(U(m)/U(m - 1) x U(l)) x fU(n)/U(n - 1) x U(l)) . 

The contributions of all these saddles to the large-w asymptotics of S^w), or 
the two-level correlation function R^u)), have been worked out in Refs. ||, |6[ 0- 
It was found that only a small number of saddle-point manifolds (p = 0,1 for 
(3 = 1, 2 and, in addition, p = 2 for f3 = 4) survive in the limit m = n = 0. 
The mechanism for termination is again the dependence on m and n of the 
(causal) analytic continuation of the volume of these manifolds. The asymptotic 
expressions for 5*2 (w) turn out to be 

(3 = 1: 1-2uj- 2 +cj- i e- 2luJ 
13 = 2: 1 - uj- 2 +uj- 2 e~ 2luJ 

(3 = 4: 1 - \u~ 2 + (tt/2) crV 3 *'' + (2w)- 4 e" 4 ^ . 

The real parts of these expressions agree with the asymptotic limits that are 
known, for i?2(w), from orthogonal polynomials |llj or the supersymmetric me- 
thod p|. It is reasonable to expect that this agreement is not accidental but 
extends to all orders in the asymptotic expansion in In other words, we 

conjecture that the fcrmionic replica trick in its recent elaboration is perturba- 
tively equivalent to the supersymmetric method. Also, note that the expression 
for (3 = 2 is not just asymptotic but is exactl We will elaborate on this later. 

What about the bosonic replica trick? For the one-point function we saw that 
there is a basic asymmetry between bosonic and fermionic replicas. The same 
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is true here. While the use of fermionic replicas for Z m ^ n {u) leads to a compact 
symmetric space X mtfl , bosonic replicas for Z_ TO ,_ n (w) lead to a noncompact 
analog space, Y m/n (which is still a Riemannian symmetric space, albeit of non- 
compact type.) The pairing of compact spaces with their noncompact analogs 
is given in the following list: 

(3=1: Sp(2?7i + 2n)/Sp(2m) x Sp(2n) 0(2m, 2n)/0(2m) x 0(2n) 

(3 = 2: U(m + n) /U(m) x U(n) U(m, n) /U(m) x U(n) 

/? = 4 : 0(2m + 2n) /0(2m) x 0(2n) Sp(2m, 2n) /Sp(2m) x Sp(2n) . 

The integral representation for Z_ m ._„(o;) in the universal low-frequency limit 
is still of the form 

Z- m ,-n(w) = y exp ^-^TrS 3 Qj . 

Moreover, the integration manifold F m .„ is still generated by acting with the 
noncompact large group G on the diagonal matrix S3 by conjugation: Q = 
gYi^g -1 {g g G), and we still have the invariant saddle point Qo = S3. The 
large-w asymptotic expansion around this saddle point of Y m>n is known JUsf to 
give the same results as the expansion around Qo = S3 for X m>n , at to = n = 0. 
What is different is that Qo now is the only saddle point. There are no others 
that lie directly on Y m>n , although there do exist saddle-point manifolds that can 
be reached from Y m>n by moving a finite distance along an imaginary direction 
of the real manifold Y m>n . For example, for (3 — 2 we set A = Ei, m+ i +^ m +i,i € 
LieU(m, n) - where E p>q is the matrix whose entries are zero everywhere except 
on the intersection of the p-th row with the g-th column where the entry is unity 
- and follow the "imaginary" orbit Q(it) = e lt T,^e~ ltA . At t = tt/2 we reach 
Q(i7r/2) = Q\. By applying the symmetry group H = U(m) x U(n) to Qi, we 
arrive at an analog of the p = 1 saddle-point manifold of the fermionic variant. 
If one evaluates the integrand on this saddle-point manifold, one gets a factor 
e — t(m+n— 2)w_ p or m _ n _ q ; j-jjjg g rows as w moves into the lower half of the 

complex plane, which is unphysical. Hence, this saddle-point manifold must be 
discarded, and we are back to the statement jL4| that the large-w limit of the 
function Z m ^ n (uj) for negative to and n is exhausted by a single saddle point 
(p = 0). In summary, while the fermionic replica trick suffers from an excess of 
saddle points, the bosonic version has "too few" . 

Although the recipe of selecting good saddle points in the fermionic for- 
mulation appears to work, it does need further justification. A more proper 
procedure is to combine the information from bosonic and fermionic replicas 
and write down a bona fide analytic continuation for Z m ^ n {iS) away from the in- 
tegers. The proposal made in |7) is again to extend the sum over saddle-point 
manifolds p all the way up to infinity. Doing so, however, we run into the same 
difficulty we analysed in considerable detail for the one-point function. We shall 
not repeat this analysis here, but only summarize the facts: for noninteger to 
and/or n, the coefficients in the sum over p grow in a factorial manner, so that 
the sum diverges and does not define an analytic continuation of Z m>n {ui) to 
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m, n ^ Z. In the absence of a well-defined analytic continuation we, of course, 
have no mathematical control on the limit — \\m. m ^ n ^Q(mn)~ 1 d 2 Z m ^ n / dw 2 . 

My interpretation of the divergence is the same as before. For a large but 
fixed value of u, the saddle-point approximation to Z m . n (u>) eventually breaks 
down (or at least it does for [3 = 1,4) when m or n becomes too large. Indeed, 
for large m = n, we should introduce standard polar coordinates on X m , n to 
cast Z m ^ n {u>) in the form of a Coulomb gas partition function, and look for a 
mean field of the Coulomb gas. Analytic continuation in the particle number 
m = n of the Coulomb gas looks benign. We therefore expect the divergent 
sum over saddle-point manifolds p to get cut off at large values of m, n. The 
existence of such a cutoff implies that there exists a nonperturbative correction 
which we are missing when writing the sum over saddle-point manifolds p. For 
large u, this correction will not affect the extrapolation to m = n — 0. On the 
other hand, for small u> the correction may cause an uncontrollable error. 

The remaining question is: what happens beyond perturbation theory? Is 
it reasonable to expect that replicas will penetrate the nonperturbative small-a; 
regime? This is hard to answer, but one warning can be issued with certainty. 
A series expansion obtained by saddle-point approximation is almost always 
asymptotic, which means the series diverges. What happens is that for a fixed 
value of the expansion parameter, 1/lu in our case, the series initially becomes 
a better approximation with increasing order, but eventually turns away and 
explodes. Without knowing the exact answer, it is often difficult to locate 
the turning point. For this reason, asymptotic expansions hardly deserve to 
be called "systematic". The only safe way of using an asymptotic series is to 
fix the order of approximation and then lower the expansion parameter 1/lu 
accordingly, so as to make the error term negligible. This limits the usefulness 
of asymptotic series in practice. To do better, one needs to establish Borel 
summability and carry out Borel resummation. This is already nontrivial for 
integrals with one saddle point, and becomes much more difficult when two 
or more saddle points are involved. On the other hand, not all saddle-point 
approximations are asymptotic, as is demonstrated by the case of the two-point 
function for (3 = 2. We will explain in the next section what is special about 
that example. 

The nonperturbative ambiguities of the replica trick for disordered electron 
systems are not restricted to the zero-dimensional limit but also affect the replica 
field theory for d-dimensional systems. Let me finish the section with an example 
demonstrating this. 

Disordered two-dimensional electrons in a strong magnetic field display the 
(integer) quantum Hall effect. A replica field theory for such systems was pro- 
posed long ago by Pruisken |l(| . The theory is a compact nonlinear a model 
with Lagrangian 

L = ^Tr d^Q d„Q + ^e^Tr Q d^Q d v Q , 
o o 

where the field Q = £/£ 3 [/ _1 parametrizes U(m + n)/U(m) x \]{n). Without 
loss in the zero replica limit, we may take m = n. The coupling constants 
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o xx and a xy are interpreted as the conductivities (dissipative and Hall) of the 
nonintcracting electron gas. The Hall conductivity a xy determines the current 
response transverse to an applied electric field. In particular, it determines the 
orientation of the current flow around the boundary of a finite sample |l7| . 

Now comes a disaster. Use of the isomorphism Q i— > — Q in Pruisken's 
Lagrangian changes the sign of the term e^Tr Qd^Qd u Q. Thus the field theory 
coupling o~ xy is defined only up to a sign. For m = n = 1, where the above field 
theory (with a xy = S) is known fl8|| to describe antiferromagnetic quantum spin- 
s' chains, this indeterminacy makes sense. Antiferromagnets carry no sense of 
orientation. But for n — 0, o~ xy is supposed to be the physical Hall conductivity, 
so the field theory makes the nonsensical prediction that the orientation of the 
current response of a quantum Hall sample is ill-determined?! 

Of course, this is just another manifestation of the causality problem we have 
been emphasizing all along. Can it be cured? Pruisken's field theory is known 
to flow to the strong coupling (or small a xx ) regime, where the dominant field 
configurations form a dense gas of interacting instantons and anti-instantons. 
In such a nonperturbative soup of saddle points and almost-saddle points, it 
seems quite hopeless to try and disentangle the good configurations from the 
bad ones. 



4 Semiclassical exactness for (3 = 2 

We have seen that the fermionic replica trick, augmented by an inspired recipe 
for selecting good saddle points, for (3 — 2 already gives the correct answer 
for the universal limit of S^w) when the saddle-point approximation is carried 
out to leading order in the small parameter 1/uj. Apparently, the saddle-point 
approximation in this case is exact, not just approximate! We call integrals 
where this miracle happens "semiclassically exact" . In the present section we 
will explain the mathematical basis underlying the phenomenon of semiclassical 
exactness: the Duistermaat-Heckman theorem. This is a celebrated result in 
symplectic geometry JH], [H| ^2| and is included here for the convenience of 
the targeted reader. A few references congenial to physicists are [|3 ( 24. p5| . 



Let M be a symplectic manifold of dimension 2/. In physics such a manifold 
is called a phase space with / degrees of freedom. Simple examples are the two- 
sphere S 2 or the two-torus T 2 . The example we are particularly concerned with 
is the coset space U(m + n)/U(m) x U(n), which has dimension 2/ = 2mn. 
A symplectic manifold M comes with a symplectic structure, i.e., a closed and 
nondegenerate two-form (I. For the two-sphere, (I is the solid angle, expressed in 
spherical polar coordinates 0, <j> by sin d9 Ad(j>. For the space U(m + n) /U(m) x 
U(n) parametrized by Q, the symplectic structure is (8i) _1 TrQdQ A dQ. By 
Darboux's theorem, there exist local coordinate systems made up of positions 
q l and momenta pi, such that (I takes the canonical form (I — J2i=i ^Pi A dq l . 
Let fi^ = Ilf=i dpi A dq 1 denote the Liouville form on M. In the case of interest, 
(It agrees with the volume element dQ up to a multiplicative constant. 

The Duistermaat-Heckman theorem in the formulation given below makes a 
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statement about integrals of the form 



Z(u)= [ n f e lujH 

JM 



where the function H : M — > M is required to possess the following property. 
Viewing H as a Hamiltonian function, we have the canonical equations of motion 
of Hamilton mechanics: 

i dH dH 
dpi dq l 

The solutions of these equations are conveniently assembled into a map MxK-> 
M, (x,t) i— > ipt(x), called the phase flow of H. The crucial property we shall 
require of H is that there exist some Riemannian metric on M which is preserved 
by the flow ipt- 

Before stating the Duistermaat-Heckman theorem, let us look at a few ex- 
amples where this condition is satisfied. The simplest one is provided by the 



Hamiltonian for a classical spin in a magnetic field 23 . The phase space M in 
this case is the two-sphere S 2 , and the motion is precession of the spin around 
the axis of the magnetic field. The frequency of precession, the Larmor fre- 
quency, does not depend on the tilt angle between the magnetic field axis and 
the axis of the spin. Thus the phase flow tpt is simply uniform rotation around 
the field axis, with the rotation angle being a linear function of time. Clearly, 
this flow preserves the natural metric d9 2 + sin 2 9 dtp 2 of S 2 . 

The classical spin in a magnetic field is but one of a large family of integrable 
Hamiltonian systems with the same property. A brief sketch is as follows. Let G 
be a compact semisimple Lie group with maximal abelian subgroup (or maximal 
torus) T. We fix some regular element A of the Lie algebra Lie(T), and consider 
the adjoint orbit Ad(G)A consisting of elements UAU^ 1 . The adjoint orbit 
comes with a natural Hermitian structure, i.e., a symplectic structure f2 as well 
as a metric tensor g, which we refrain from writing down here as this would 



lead us just a little too far (for these details consult Appendix 2 of 26 1). Note 
that we get the phase space S 2 of the classical spin by putting G — SU(2) and 
A = 03. Now we pick a quadratic form, say (X, Y) = Tr XY, on the Lie algebra 
Lie(G), and we fix another regular element B of Lie(T) to form the function 
H = Tr BUAU^ 1 . We may view the adjoint orbit Ad(G)^4 as a phase space, 
with H being a Hamiltonian function on it. The phase flow of this Hamiltonian 
turns out to be ipt : UAU -1 i— > e tB UAU~ 1 e~ tB , and this flow preserves the 
Riemannian metric g of the adjoint orbit. We are mentioning this class of 
examples because it will be seen to give rise to the Itzykson-Zuber formula. 
Our example of interest is the function H — Tr/iQ/2, with 

h = diag(h + , h-) = diag(/i+i, h +m ; h-i, /i_„) 

a diagonal real matrix, on the symplectic manifold U(m+n) /U(m) x U(n) . (This 
example can actually be regarded as a degenerate limit of the above general 
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family.) To construct its phase flow in an elementary fashion, we express Q as 



/ 1 -2BB^_ 2BVT~^WB\ 
Q ~ \2BWl - BB^ -l + 2B^B ) 



where B is a complex m x n matrix with range < B^B < 1. The symplectic 
structure and the Hamiltonian for this choice of coordinates take the form 

n = (8i)" 1 TrQdQAdQ = r 1 TrdSAdS t , 

H = TrhQ/2 = const - Tr h+BB f + Tr h-B^B . 

The phase flow if) t is then readily seen to be B i— > Q lth + Be~ %th ~ , or in the original 
representation, 

ipt ■ Q <-> e lth Qe~ uh . 

Clearly, this preserves the natural Riemannian metric TrdQ 2 . 

We now come to the Duistermaat-Heckman theorem. The theorem is rooted 
in equivariant cohomology, and its best version relies on no more than Stokes' 
formula for differential forms that are equivariantly closed for a compact group 
action. In particular, no reference to any metric structure is needed. Neverthe- 
less, to avoid mathematical overload, we shall permit ourselves the luxury of 
using a metric. 

Theorem. Let (M, Vl, H) be a Hamiltonian system with compact 2/-dimen- 
sional phase space M, and let its phase flow preserve a Riemannian metric on 
M. Then the integral J M We 1,11 localizes on the critical set of H . 

Remark. What the theorem is saying is that the integral is completely 
determined by the values of the function H, and a finite number of derivatives 
thereof, on the set of solutions of dH = (the saddle points or saddle-point 
manifolds). In other words, the stationary-phase approximation is exact. 

Proof. We will employ a method of proof originally due to Bismut |2?j (see 
also Refs. |2^, ^]), which uncovers and takes advantage of a hidden supersym- 
metry. Let x 1 (i — 1, ...,2/) be a system of local coordinates of M, in which 
the symplectic structure is expressed by ft = ^ilij(x) dx l A dx^ (summation 
convention). We supersymmetrically extend the phase space M by introducing 
for every x l a corresponding anticommuting coordinate £ 4 . Consider then the 
expression 

where / means integration with the flat Berezin form da: 1 . . . dx 2 f d^i . . . dgsf. By 
using canonical coordinates, it is not difficult to see that the Fermi integral, i.e., 
differentiation with respect to the anticommuting variables and multiplication 
by dec 1 . . . dx 2 ? , produces the Liouville form . Hence we have 



J exp ( + = J fl f e iH 
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Next, using the closedness of f2 ^dCljk/dx 1 = 0), one verifies that the 
exponent 5 = —kSljk£'£ k + iH has the property of being annihilated by 

r-, ^ d ,^ ik . dH d 

£> = . -)- ^ri- 3 - 

dxi dxi d£ k 

D is a first-order differential operator of odd type, and is analogous to the 
BRST operator that plays a central role in the functional integral quantization 
of nonabelian gauge theories ]3(|. The key step now is to deform the integral 
by a parameter t: 

' nfe* H = /e s = [ e s+tDX 



M 

2 • 



where the function A is constrained by D A = and will be specified presently. 
The Berezin integral on the right-hand side has the crucial property of being 
independent of t. Indeed, by Taylor expanding, 



Je s+tDX = J c s (l + tD\+\t 2 {D\) 2 + 



and using D 2 X — and partial integration in conjunction with DS = 0, one sees 
that the terms of linear and higher order in t all vanish. To reap full benefit from 
the i-independence of the integral, we pick a well-chosen Ricmannian metric g 
and set 

It can be shown that this expression satisfies the condition D 2 X — if and only 
if the metric g is invariant under the phase flow of H . By assumption, such a 
metric g exists. The last step is to take the parameter t to infinity. Since the 
number part of 

fjflfifT 

DX = ^ dxldx^ + ° m = - ldHl2 + ° m 

is negative definite, this limit localizes the integral J" ^e lH = J e s+tDA onto 
the critical set of H, i.e., the solutions of the saddle-point equation dH = 0. 
This concludes the proof. 

If the saddle points of H are isolated, one can easily push the calculation 
further and write down an explicit formula for the integral as a discrete sum: 

/ n f e lH = (2vr) / V ^sgnCHess.g) QlH[x) 
Jm x:d f^ =0 ^\Bct(Ress x H)\ ' 

where the phase of the contribution from each saddle point is determined by the 
signature of the Hessian of H at x. For the classical spin in a magnetic field, 
this gives 

e dn = 2n — ijlj — ' 

n 2 = l 
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the correctness of which is easily verified by direct computation. Less trivial inte- 
gral formulas result on taking for the Hamiltonian system (M, fi, H) the adjoint 
orbit of a compact Lie group G, with the Hamiltonian being H = Tr AUBU -1 . 
If both A and B are regular elements in a Cartan subalgebra of G, one obtains 

Mm 



j e ^AUBU-- du=ij){A)p{B) yl (-!)I"V V 



where p(A) — ]~[r»->n is a product over positive roots, the sum runs over the 
Weyl group of G, and \w\ denotes the parity of w. (Regularity of A means 
that p(A) is nonzero.) dU is a suitably normalized Haar measure of G. Note that 
the integral can and should actually be understood as being over the symplcctic 
quotient G/T. That the factor of i has disappeared from the exponent is of no 
concern. Looking back at the proof of the Duistermaat-Heckmann theorem, we 
see that i didn't play any role: the underlying principle is localization on the 
critical set, not stationary phase. 

The above formula is due to Harish-Chandra (1957) who proved it by a 
quite different method, namely by computing the radial parts of G-invariant 
differential operators pl}| . Specializing to the case G — U(iV) (or SU(iV), it 
doesn't matter) and setting A — diag(A l7 . . . , A N ), B = diag(f?i, . . . , B N ) we 
get 

^aubu- 1 du _ Det(e iBj )i, j= i,...N 



/u(w) Ui^-AjKBi-Bj)' 

which is known in physics as the Itzykson-Zuber formula. 

Wc finally understand the semiclassical exactness of the replica integral for 
the (3 = 2 two-point function in the universal limit. The nonlinear tr model man- 
ifold U(to + ti)/U(to) x U(n) is a 2mn-dimensional phase with symplectic struc- 
ture (8i)" 1 TrQd(5AdQ, and the flow of the "Hamiltonian" H = wTrS 3 Q/2 pre- 
serves a Riemannian metric g = TrdQ 2 . Therefore, the Duistermaat-Heckman 
localization principle applies, and the integral J exp(iwTr S 3 Q/2) dQ is evalu- 
ated exactly by saddle-point approximation. 



5 Conclusion 

The replica trick for spectral correlations of disordered electron systems is math- 
ematically ill-founded because analytic continuation of the generating function 

Z m>n (u, v) = (Det m {u - H) V>ct n {v - H)} 

away from the integers m,n is not unique. In benign cases, such as the famous 
Selberg integral, uniqueness is guaranteed by a boundedness property as spelled 
out by Carlson's theorem or Carleman's theorem ]32] ). It appears that no such 
property is available in the present case. (However, it has been suggested to me 
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[ p3| that this problem might be overcome by considering a generating function 
of the positive real form 

(Det m ((^ 1 - Hf + e 2 ) Det n ((E 2 - H) 2 + e 2 )) .) 

For positive m, n (the fermionic replica trick) the generating function is 
analytic in u and v, which means that the entire information about causal- 
ity of Green's functions is missing. (Of course, the enlightened user knows 
about causality, but fermionic replicas by themselves do not.) In recent work by 
Kamenev and Mezard, and by Yurkevich and Lerner, causality was introduced 
into the formalism by selecting a well-chosen set of saddle points of the compact 
nonlinear a model based on fermionic replicas. In this way, the known large- 
frequency asymptotics of the level correlation functions for disordered electron 
systems was reproduced for all cases (3 = 1,2,4. This motivates the conjecture 
that the fermionic replica trick (augmented by the selection of "causal" saddle 
points) is perturbatively equivalent to the supersymmetric method. 

Close inspection of Refs. || ^, Q reveals that the procedure used there is 
mathematically uncontrolled, as it involves a nonexistent analytic continuation. 
The fact of the matter is that the sum over saddle-point manifolds p = 0, 1, ... oo 
diverges for noninteger m, n. We have argued that this divergence is likely to be 
cut off by a nonperturbative mechanism. Although such a cutoff is immaterial 
for the analytic continuation to m = n = for large frequency uj, it may 
introduce uncontrollable errors for small ui. In any case, the method as it stands 
is limited to the perturbative regime, as it relies on stationary-phase evaluation 
of integrals, which requires a small parameter 1/uj. 

A miraculous exception is Wigner-Dyson statistics for (3 = 2, where a hidden 
supersymmetry, namely equivariant cohomology and the localization principle 
underlying the Duistermaat-Hcckman theorem, localizes the nonlinear a model 
integrals on the saddle points. Put differently, the leading-order stationary- 
phase approximation in this case is exact, regardless of whether uj is large or 
small. The same mechanism is at work in the supersymmetric formulation. 
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